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CLAIMS 

1. A method of determining hysteresis of a process device in a 
process, the method comprising 

collecting sample pairs of signals representing a control input and a 
5 process output of a process device, 

calculating the hysteresis of the process device on the basis of the 

collected sample pairs, 

characterized by 

each collected sample pair comprising an average control input and 
10 an average output for a predetermined collection period, 

screening (53) the sample pairs suitable for further processing at 
least according to the magnitude and direction of the relative change of the 
average control input of each sample pair, 

grouping (53) the sample pairs Into a first and a second group at 
15 least according to the direction of the relative change of the average control 
input of each sample pair, 

calculating (54) a first and a second control/output characteristic 
curve by means of the screened sample pairs of the first and the second 
group, respectively, 

20 determining the hysteresis of the process device as a distance be- 

tween the first and the second characteristic curve at least at one point. 

2. A method according to claim 1, characterized in that the 
signal representing the process output comprises one of the following: a meas- 
ured process variable (y), or a set point (r) of a process variable in process 

25 control. 

3. A method according to claim 1 or 2. characterized in that 

the screening (53) comprises 

calculating the change of average control input for each sample pair 
in relation to the average control input of the previous sample pair and the di- 

30 rection of the change, 

if the calculated change is smaller and has occurred in the same di- 
rection as the change calculated for the control input of the previous sample 
pair, the sample pair is selected for further processing; othenA/ise the sample 
pair is rejected. 
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4. A method according to claim 1, 2 or 3, characterized in 
that the direction of the change is determined by calculating the sign of the 
change, in which case the sign is positive if the control input is Increasing, and 
negative if the control input is decreasing. 
5 5. A method according to claim 1, 2, 3 or4, characte rized in 

that the screening comprises 

selecting a sample pair UmO). YmO) for the calculation of the charac- 
teristic curve if the following two conditions are fulfilled: 

lAu„(j)l<KO-l)| 
signiAu^ (J)) = signihu^ (J - 1)) 

where 

UmG) "'s the average control input of sample pair j, 
YmC) is the average measured output of sample pair j, 
15 j is an integer Index, 

AUmO)=UmG)-UmO-1) 

sign(AUm) calculates the sign of AUm- 

6. A method according to claim 5, characterized in that the 
sample pairs are grouped into a first and a second group on the basis of the 

20 following rule: 

if the sign of AUm is positive, the sample pair belongs to the first 

group. 

if the sign of AUm is negative, the sample pair belongs to the second 

group. 

25 7. A method according to claim 1,2, 3 or 4, characterized in 

that the collecting comprises 

taking a momentary sample pair (ys, Us) from the control input (Um) 
and measured output (y^) at certain intervals, the interval being preferably in 
the order of one or more seconds, 

30 calculating the average sample pair (Um, ym) from the momentary 

sample pairs (ys. Us) taken during the collection period, which is preferably in 
the order of one or more minutes. 

8. A method according to any one of the preceding claims, char- 
acterized by calculating a reliability value for hysteresis by a function, 

35 which includes the following information as parameters 
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the number Ny1 of average sample pairs (Um. Ym) belonging to the 

first, ascending characteristic curve, 

the number Ny2 of average sample pairs (Um. Ym) belonging to the 

second, descending characteristic curve, 
5 the number Noy1 of sample pairs which belong to the first charac- 

teristic cun/e but are below the second characteristic curve, and 

the number Noy2 of sample pairs which belong to the second char- 
acteristic curve but are above the first control curve. 

9. A method according to claim 8, c h a r a c t e r i z e d in that 

0 the reliability value L is calculated by function L=max(0.1 -Noy1 /Ny1 - 

Noy2/Ny2). in which case L may obtain values from 0 to 1 , where L=1 is com- 
pletely reliable and L=0 is completely unreliable. 

10. A method according to any one of the preceding claims, 
c h a r a ct e r i z e d in that the calculation of each characteristic curve com- 

15 prises 

dividing the control area of the control input into bins uo(1)...Uo(nbin). 
where Uo(1)...Uo(nbin) is the location of the bins on the u-axis representing the 
control Input and bin is the number of bins, in which case the values 
Yo(1)...yo(nbin) included in the bins represent output values on the y-axis de- 
20 scribing the output, and value pairs uo(1 )/yo(1 ) Uo(nbin)/yo(nbin) define a char- 
acteristic curve, 

updating the values yo(1)-..yo(nb.n) of the bins uo(1)...Uo(nbin) by 
means of the average sample pairs using a predetemiined weighting function. 
11. A method according to claim 10, c h a r acte ri z ed by 
25 updating two bins bn. bn-i according to each sample pair, the bins 

being selected so that the following condition is fulfilled 



fc„< (n,,„-l) + l<6„.. 



where 

Um is the minute mean value of the control of sample pair j, 
3Q is the minute mean value of the measurement of sample pair j, 

n = 1...bin, 

Umin and Umax are the minimum and the maximum of the control area. 



respectively, 

updating the values yo(bn) and yo(bn+i) of the selected bins bn. and 
35 bn-1 as follows 
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10 



15 



25 



ncKh)yM + w,y„ 
^"^"^ nct(b„) + w, 

ncf(6,.„)3^o(&„^i) + W2ym 

where 

nct(1)...nct(nw„) each represents the number of updates (con 
trol/measurement pairs used) of each bin, 
Wn and Wn+1 are weighting coefficients 



_^ \u^-Unin-(P„-l)Ust\ 



w^i=i 

Ust is the distance between the bins 



W - 1 



updating the numbers nct(bn) and nct(bn-1) of updates of the se- 
lected bins bn and bn-1 as follows 
nct(b„) = a*nct(bn) + wi 
nct(bn+i) = a*nct(bn+i) + W2 
20 where a is constant. 

12. A system for determining hysteresis of a process device in a 

process environment, the system comprising 

means for collecting sample pairs of signals representing a control 

input and a process output of the process device, 

0 h a r a c t e r i z e d in that each collected sample pair comprises 
an average control input and an average control output for a predetemiined 
collection period, and the system comprises 

means (53) for screening sample pairs that are suitable for further 
processing at least according to the magnitude and direction of the relative 
30 change of the average control input of each sample pair. 
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means (53) for grouping the sample pairs into a first and a second 
group at least according to the direction of the relative change of the average 
control input of each sample pair, 

means (54) for calculating a first and a second control/output char- 
5 acteristic curve by means of the first and the second screened sample pair (Um, 
ym), respectively, 

means (55) for detemiining the hysteresis of the process device as 
a distance between the first and the second characteristic curve at least at one 
point. 

10 13. A system according to claim 12, characterized in that 

the signal representing the process output comprises at least one of the follow- 
ing: a measured process variable (y), or a set point of a process variable (r) in 
process control. 

14. A system according to claim 12 or 13, characterized in 

1 5 that the screening means comprise 

means (53) for calculating the change of the average control input 
for each sample pair in relation to the average control input of the previous 
sample pair and the direction of the change, and for selecting a sample pair for 
further processing if the calculated change is smaller than the change calcu- 

20 lated for the conti-ol input of the previous sample pair and has occured in tiie 
same direction; othenwise the sample pair is rejected. 

15. A system acconding to any one of claims 12 to 24, c h a r a c - 
t e r I z e d in that the system comprises means (55) for calculating a reliability 
value for hysteresis by a function which includes the following infonnation as 

25 parameters 

the number Nyl of average sample pairs (Um. ym) belonging to tiie 

first, ascending characteristic curve, 

the number Ny2 of average sample pairs (Um, ym) belonging to the 
second, descending characteristic curve, 
30 the number Noyi of sample pairs which belong to the first charac- 

teristic curve but are below the second characteristic curve, and 

the number Noy2 of sample pairs which belong to the second char- 
acteristic curve but are above tiie first control curve. 

1 6. A system according to claim 14, characterized in that 
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the reliability value L is calculated by function L=max(0.1-Noy1/Ny1- 
Noy2/Ny2), in which case L may obtain values from 0 to 1. where L=1 is com- 
pletely reliable and L=0 completely unreliable. 

1 7. A program product which includes a program code, which im- 
5 plements the method steps of a method according to any one of claims 1 to 1 0 
when it is run on a computer or the like. 



